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Рассматривается возможность применения линейных нейронных сетей при решении задач прогноза технологических показате- 
лей нефтедобычи. Исследуются эффективность и помехоустойчивость нейросетевых алгоритмов. 


В опубликованных ранее работах [1-4] в основ- 
ном рассматривались методологические проблемы 
решения задач нефтепромысловой геологии с ис- 
пользованием нейросетевых алгоритмов. При этом 
выбор архитектуры сети осуществлялся исходя из 
принципа ее универсальности, а основные расчеты 
осуществлялись с использованием многослойной 
нейронной сети обучающейся по методу "Ьаск ргор- 
а§аііоп" и способной решать нелинейные задачи. 
Однако, на наш взгляд, не следует пренебрегать и 
более простыми линейными нейронными сетями. 

Зачастую можно встретить такую ситуацию, ког- 
да задача на первый взгляд кажущаяся сложной и не- 
линейной, на самом деле может быть успешно реше- 
на линейными методами. Таким примером может 
служить задача прогноза показателей нефтедобычи 
по данным ежемесячных регламентных наблюдений. 

Цель работы - исследование точностных харак- 
теристик и помехоустойчивости нейросетевых ме- 
тода прогноза показателей нефтедобычи с исполь- 
зованием линейных нейронных сетей. 

На языке нейронных сетей линейная модель 
представляется сетью без промежуточных слоев, 
которая в выходном слое содержит только линей- 
ные элементы (т.е. элементы с линейной функцией 
активации). Веса соответствуют элементам матри- 
цы, а пороги - компонентам вектора смещения. Во 
время работы сеть фактически умножает вектор 
входов на матрицу весов, а затем к полученному 
вектору прибавляет вектор смещения. Одним из 
наиболее распространенных методов обучения ли- 
нейной нейронной сети является стандартный ал- 
горитм линейной оптимизации, основанный на 
псевдообратных матрицах, а также алгоритм с 
адаптивной настройкой шага обучения [5, 6]. 


Исследования проводились с использованием 
методов имитационного моделирования по схеме 
замкнутого численного эксперимента. Известно 
[8, 9], что все разнообразие моделей можно условно 
свести к трем типам: физическая (натурная) мо- 
дель, аналоговая и математическая. Примером ма- 
тематической модели в нашем случае могут служить 
системы тех или иных уравнений, описывающих 
процессы фильтрации жидкостей в пористых сре- 
дах [7, 9]. Однако, учитывая то, что основной целью 
исследований являются не сами процессы фильтра- 
ции, а нейросетевые методы, можно ограничиться 
достаточно простой физической моделью, описы- 
вающей процессы перераспределения давления в 
пласте под воздействием работы системы добываю- 
щих и нагнетательных скважин, а именно на моде- 
ли одномерной однофазной фильтрации жидкости. 

В основе данной модели лежит широко применяе- 
мая формула Тэйса [9], согласно которой понижение 
давления в А р(гф) в любой момент времени 1 в точке 
пласта, расположенной на расстоянии г от возмущаю- 
щей скважины определяется до и после остановки сква- 
жины в момент времени Т следующими формулами: 
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АР{Г ' 1) = ІЙё (ЕІ( 4^ +ЕІ( 4^Г) }) ПРИ ' > Г ’ П 

где Еі(-х) = | Р^-сІи - интегральная показательная 


функция. 

В данной формуле 0 - среднесуточный дебит 
скважины, /и - вязкость жидкости, к - проницае- 
мость пласта, Ъ - толщина пласта, х ~ коэффици- 
ент пьезопроводности. 
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Формула (*) позволяет легко рассчитать поле 
давлений в любой точке пласта, создаваемое рабо- 
тающей скважиной. В свою очередь, для системы 
работающих в различных режимах скважин, давле- 
ние в любой точке пласта может быть рассчитано 
как сумма полей давлений от отдельных скважин. 

Данная расчетная схема была реализована на 
ПЭВМ. При ее численной реализации учитывался 
тот факт, что по результатам расчетов должна быть 
построена обучающая выборка, по составу совпа- 
дающая со стандартным набором параметров, со- 
держащихся в соответствующих базах данных за- 
полняемых в процессе разработки реальных место- 
рождений. 

Сделаем ряд замечаний относительно методов 
обучения нейронных сетей. Как известно, для обу- 
чения нейронной сети требуется обучающая вы- 
борка, состоящая из отдельных примеров, каждый 
из которых представляет собой задачу одного и то- 
го же типа с соответствующим набором входных 
параметров и ответом. 

В нашем случае моделируемая согласно (*) сис- 
тема полностью определяется координатами добы- 
вающих и нагнетательных скважин, режимами их 
эксплуатации, временем их работы и простоя, гео- 
логическими характеристиками пласта. Отметим, 


что конкретный вид обучающей выборки зависит 
от типа решаемой задачи, но согласно [1^1] их 
можно условно разделить на два типа - вертикаль- 
ные и горизонтальные (табл. 1,2). 

В табл. 1 приняты следующие обозначения: 
г) (/=1,2, ... ,т ) -время отбора информации, м-коли- 
чество примеров в обучающей выборке для всех сква- 
жин, Іп)ц (і=\,2,...,т к , - объемы закачки 

для нагнетательных скважин, Л) - количество нагне- 
тательных скважин, Ті и (і=\,2,...,т к ,і=1,2,...,Щ -вре- 
мя работы нагнетательных скважин, Тр ч (і=\,2,...,т к , 
]=\,2,...,Щ - время работы добывающих скважин, 
Ы р - количество добывающих скважин. 

Отметим, что в общем случае в поля обучения 
могут быть включены данные о пластовых и забой- 
ных давлениях, глубинах погружения насосов, ди- 
намическом уровне, технологических характерис- 
тиках оборудования, геологических характеристи- 
ках пласта и т.д. 

Согласно [1^1] горизонтальные выборки в ос- 
новном используются при решении задач исследо- 
вания интерференции скважин и разделения эф- 
фектов от проводимых геолого-технических ме- 
роприятий, а также для решения задач прогноза 
суммарной добычи нефти для выделенных фраг- 
ментов (ячеек) месторождения. 


Таблица 1. Структура простейшей горизонтальной обучающей выборки 
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Таблица 2. Структура простейшей вертикальной обучающей выборки 
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Их существенный и неустранимый недостаток 
- количество полей обучения возрастает прямо 
пропорционально количеству скважин и, следова- 
тельно, область их применения ограничена иссле- 
дованиями выделенных ячеек месторождения с ко- 
личеством скважин порядка 30-50, в то время как 
число скважин в отдельных месторождениях может 
быть порядка нескольких сотен и даже тысяч. 

Второй существенный недостаток подобных 
выборок - высокие требования к входной инфор- 
мации, т.к. наличие пробелов в данных по отдель- 
ным скважинам для различных моментов времени 
приводит к тому, что соответствующие строки в 
обучающей выборке приходится отбрасывать, либо 
заменять их приближенными значениями. 

В свою очередь вертикальные обучающие вы- 
борки предназначены для решения задач прогноза 
добычи нефти по отдельным скважинам и решения 
задач классификации [1-4]. Они менее чувствитель- 
ны к наличию пробелов в базах данных, но абсолют- 
но не приспособлены к решению таких задач, как 
исследование интерференции скважин. В табл. 2 
приведена структура простейшей обучающей вы- 
борки. В ней приняты следующие обозначения: 
(і (/=1,2,... ,т) - время отбора информации т - коли- 
чество примеров в обучающей выборке для всех 
скважин, ІЩф (/=1,2,. -приемистос- 
ти либо объемы закачки для нагнетательных сква- 
жин, Ы, - количество нагнетательных скважин, 
Тіф (/=1,2,. ,.,т ь у=1,2 ,. ,.,Щ - время работы нагнета- 
тельных скважин, Тр ф (/=1,2, ...,пг к , у- 1 , 2 - 
время работы добывающих скважин, Ы р - количество 
добывающих скважин, 0 Ч , /=1,2,...,/я ь - 

дебиты добывающих скважин, X, (і=1,2,...,Щ - ко- 
ординаты, либо номера добывающих скважин. 

Построим простейшую горизонтальную выбор- 
ку, где полями обучения нейронной сети будут: те- 
кущее время, координаты, время работы и объем 
закачки для нагнетательных скважин, время рабо- 
ты добывающих скважин и их координаты, а в ка- 
честве прогнозируемого параметра выберем сум- 
марный среднесуточный дебит добывающих сква- 
жин, либо дебит для одной из добывающих сква- 
жин. Очевидно, что при такой постановке задачи 
мы сознательно ограничиваем набор входных дан- 
ных (количество полей обучения). Более того, и 
при выборе периода обучения сети, мы также огра- 
ничивали период обучения до 2,5... 3 лет. 

Причиной этого служит тот факт, что в реаль- 
ной ситуации в процессе разработки нефтяного 
пласта его коллекторские свойства могут изменять- 
ся в весьма значительной степени. Следовательно, 
выбрав большой временной период обучения, мы 
столкнемся с тем фактом, что примеры, содержа- 
щиеся в обучающей выборке, не будут однородны- 
ми по содержанию. 

Иными словами, обучив нейронную сеть и сведя 
ошибки обучения к минимуму, мы, тем не менее, не 
сможем гарантировать того, что ошибки обобщения 
на период прогноза будут небольшими. Таким обра- 


зом, время обучения сети надо выбирать оптималь- 
ным образом. Согласно [10] этот период не должен 
превышать 2,5...3 года. Именно поэтому при расчетах 
обучающих выборок с использованием соотношений 
(*) выбирался четырехлетний расчетный период, а 
именно 3 года для обучения и 1 год для прогноза. 

Рассмотрим схему проведения замкнутого чис- 
ленного эксперимента, предназначенного для ис- 
следования эффективности нейросетевых методов 
прогноза показателей нефтедобычи. Согласно об- 
щепринятой методике [1-4, 7-8] она состоит из 
нескольких этапов, а именно: 

1. Задаются схема расположения скважин, время 
разработки и параметры пласта, характеристи- 
ки скважин. 

2. На основе соотношений (1) строится горизон- 
тальная обучающая выборка (в данном конк- 
ретном случае 36 примеров для обучения и 12 
для прогноза). 

3. Моделируется погрешность эксперимента. При 
этом в соответствующие поля обучающей вы- 
борки в период обучения вносится заданная 
случайная относительная погрешность опреде- 
ленного уровня. 

4. Строится нейронная сеть определенной архи- 
тектуры (в данном случае линейная) на ее вход 
подаются зашумленные данные, и проводится 
ее обучение. 

5. На вход обученной нейронной сети подаются 
данные за весь исследуемый период (обучение- 
прогноз) и на выходе получаем искомое решение. 

6. Рассчитываются ошибки обучения и прогноза. 
Обратимся к полученным результатам, но 

прежде отметим, что оценка точности обучения и 
прогноза рассчитывалась для четырех уровней вно- 
симой относительной погрешности, которые зада- 
вали следующим образом 0, 10, 20, 30 % для деби- 
тов и приемистостей и, соответственно, 0, 2, 4, 6 % 
для времени работы скважин. 

Цель данных расчетов - оценить помехоустой- 
чивость линейной сети в зависимости от уровня 
вносимой погрешности и количества добывающих 
скважин. 



Рис. 1 . Прогноз суммарного дебита для системы из 2-х наг- 
нетательных и 3-х добывающих скважин. Довери- 
тельный интервал 25 %. Расстояние между скважи- 
нами 400 м. Г оризонтальные выборки: 1) точные зна- 
чения суммарного дебита, 2, 3, 4) прогнозные значе- 
ния, полученные при 10, 20, 30 % уровне вносимых 
погрешностей 
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В качестве примера на рис. 1 приведем резуль- 
таты прогноза суммарного дебита. Здесь и далее 
расчеты проводились для следующей схемы распо- 
ложения скважин: в ряду скважин между двумя 
нагнетательными скважинами расположены три 
добывающих. При проведении расчетов изменя- 
лось лишь количество скважин, расстояние между 
ними и характеристики пласта. 

Анализируя представленные на рис. 1 графики 
легко заметить что для данной системы ошибки 
прогноза суммарного дебита не превышают 30 %. В 
свою очередь, аналогичные расчеты показали, что 
для системы состоящей из 4 нагнетательных и 9 до- 
бывающих скважин ошибки прогноза уже не пре- 
вышают 10 %. Более детальное рассмотрение поз- 
воляет заметить, что для горизонтальной обучаю- 
щей выборки увеличение количества добывающих 
скважин приводит к уменьшению ошибок прогно- 
за суммарного дебита. Этот эффект можно оценить 
количественно. С этой целью рассчитывались мак- 
симальные и средние оценки погрешностей прог- 
ноза для различных уровней вносимых погрешнос- 
тей и количества скважин в исследуемых системах. 

Соответствующие оценки погрешностей прог- 
ноза в зависимости от числа скважин и уровней 
вносимых погрешностей представлены на рис. 2 и 
рис. 3. Приведенные на них зависимости хорошо 
иллюстрируют эффект уменьшения ошибок прог- 
ноза суммарного дебита при увеличении количест- 
ва полей обучающей выборки и позволяют количе- 
ственно оценить среднюю и максимальную ошиб- 
ку прогноза в зависимости от уровня вносимых 
погрешностей и количества добывающих скважин. 
Отметим, что, как показали детальные расчеты, 
дальнейшее увеличение числа скважин не приво- 
дит к уменьшению ошибок прогноза. 



Рис. 2. Оценка максимальных погрешностей прогноза сум- 
марного дебита в зависимости от уровня вносимых 
погрешностей и количества скважин. Горизонталь- 
ные выборки: 1, 2, 3) погрешность прогноза при 10, 
20, 30 % уровне вносимых погрешностей 

Перейдем далее к исследованию точностных ха- 
рактеристик нейровычислений с использованием 
линейной сети и вертикальных обучающих выбо- 
рок. С этой целью поступим следующим образом: 
построим вертикальные выборки, используя уже 
полученные ранее расчетные данные. Такой прием 
позволит нам легче сопоставить результаты модели- 
рования для этих двух типов обучающих выборок и 
сделать соответствующие выводы и рекомендации. 



Рис. 3. Оценка средних погрешностей прогноза суммарно- 
го дебита в зависимости от уровня вносимых пог- 
решностей и количества скважин. Горизонтальные 
выборки: 1, 2, 3) погрешность прогноза при 10, 20, 
30 % уровне вносимых погрешностей 

Однако если ранее прогнозируемым парамет- 
ром являлась суммарная добыча, то в данном слу- 
чае эту характеристику можно определить ли ш ь 
как результат суммирования результатов прогноза 
для отдельных скважин и соответствующие оценки 
нужно строить для каждой отдельно взятой сква- 
жины. Поэтому при построении оценок качества 
прогноза и обучения мы остановились на следую- 
щем варианте. 

Максимальные оценки обучения и прогноза 
выбираются из всей совокупности прогнозных 
данных, т.е. результатов обучения и прогноза по 
всем скважинам. Аналогично можно рассчитать и 
средние ошибки обучения и прогноза. Эти характе- 
ристики были рассчитаны для ранее используемых 
систем. В отличие от результатов прогноза суммар- 
ного дебита получаемого с использованием гори- 
зонтальной обучающей выборки здесь в явной 
форме не наблюдается эффект уменьшения оши- 
бок обучения и прогноза в зависимости от числа 
скважин. Однако при последующем расчете сум- 
марных дебитов на основе данных прогноза по от- 
дельным скважинам эффект уменьшения ошибок 
прогноза наблюдается, но причина его появления 
скрыта в регуляризирующих свойствах нейросете - 
вых алгоритмов. 

Рассмотрим представленные в табл. 3 данные. 
Легко видеть, что в данном случае максимальные 
ошибки обучения для отдельных скважин могут 
превышать уровень вносимых погрешностей. При 
этом характер зависимости ошибок обучения и 
прогноза перестает монотонно зависеть о количе- 
ства скважин, а скорее определяется уровнем вно- 
симых погрешностей, Но, тем не менее, как и в 
случае применения горизонтальных выборок мож- 
но сказать, что максимальная ошибка прогноза 
суммарного дебита не превышает погрешностей 
определения полей обучающей выборки (приемис- 
тостей, дебитов) и убывает с увеличением числа ис- 
следуемых скважин. 

Таким образом, на основании результатов чис- 
ленного моделирования можно сделать следующие 
выводы: 
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Таблица 3 . Вертикальные выборки. Зависимость средних и максимальных ошибок прогноза суммарного дебита в зависимости 
от числа добывающих скважин при различном уровне вносимых погрешностей 


Погрешность в % 

Количество скважин 

3 

9 

15 

Вносимая 

0 

10 

20 

30 

0 

10 

20 

30 

0 

10 

20 

30 

Максимальная 
ошибка обучения 

2,5 

13,3 

35,8 

37,9 

2,8 

13,4 

26,9 

35,9 

4,8 

12,4 

24,6 

48,9 

Средняя 

ошибка обучения 

1,1 

5,6 

12,0 

15,4 

1,0 

4,7 

9,4 

18,8 

1,4 

4,2 

6,8 

13,3 

Максимальная 
ошибка прогноза 

2,2 

2,5 

13,9 

12,1 

2,2 

7,2 

5,0 

4,6 

3,3 

2,5 

3,5 

7,5 

Средняя ошибка 
прогноза 

0,9 

1,3 

8,8 

3,3 

1,1 

2,6 

3,3 

2,2 

1,7 

1,7 

1,9 

3,3 


1. Линейные нейронные сети при минимальном 
наборе параметров обучения позволяют эффек- 
тивно решать задачи прогноза суммарного де- 
бита для выделенных фрагментов (ячеек) неф- 
тяного месторождения. 

2. Использование как горизонтальных, так и вер- 
тикальных обучающих выборок позволяет за 
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счет увеличения количества полей обучения 
(количества исследуемых скважин) минимизи- 
ровать ошибки прогноза. 

3. Максимальная ошибка прогноза суммарного 
дебита не превышает погрешностей определе- 
ния полей обучающей выборки и убывает с уве- 
личением числа исследуемых скважин. 
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